November 26, 2007 



A primer on elliptic functions with applications in classical mechanics 
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The Jacobi and Weierstrass elliptic functions used to be part of the standard mathematical arsenal 
of physics students. They appear as solutions of many important problems in classical mechanics: 
the motion of a planar pendulum (Jacobi), the motion of a force- free asymmetric top (Jacobi), the 
motion of a spherical pendulum (Weierstrass), and the motion of a heavy symmetric top with one 
fixed point (Weierstrass). The problem of the planar pendulum, in fact, can be used to construct 
the general connection between the Jacobi and Weierstrass elliptic functions. The easy access to 
mathematical software by physics students suggests that they might reappear as useful tools in the 
undergraduate curriculum. 
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I. INTRODUCTION 

A long time ago, physics students were well trained in 
applications of elliptic functions in solving a great variety 
of problems in classical mechanics. For example, Whit- 
taker's Treatise on the Analytical Dynamics of Particles 
and Rigid Bodies^ appears to implicitly assume that the 
reader is fully conversant with the theory of elliptic func- 
tions. 

Elliptic functions rapidly fell out the standard physics 
curriculum over the past fifty years, however, and they 
are now only mentionned in passing in most stan- 
dard (modern) textbooks on classical mechanics. ^^^^^ The 
information on these mythical functions is now rele- 
gated to mathematics textbooks^*^ and mathematical 
handbooks^'^ with notations and conventions that are 
often contradictory or difficult to understand by physi- 
cists. The purpose of this paper is thus to (re)introduce 
the Jacobi and Weierstarss elliptic functions to a new 
generation of physics students through a series of stan- 
dard problems found in classical mechanics. 

Many problems in classical mechanics involve calculat- 
ing of one of the following integrals.^ The time integral 
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arises in solutions x{t) of one-dimensional problems asso- 
ciated with motion of a particle of mass m and constant 
total energy £^ in a time-independent potential U{x)^ 
where the initial condition xq may be chosen to corre- 
spond to a root of the turning-point equation E = U{x). 
The orbit integral 



e{s) 



± 



da 



V(2/x/n[i?-C/(a-i)-£VV2M] 



(2) 



mass /i, with constant total energy E and angular mo- 
mentum i = jir^O^ in a central potential /7(r), where the 
initial condition ^(0) = 5o is a turning point. The time 
integral 



m 



du 



^{2/h){l-u^)[E-V{u)] 



(3) 



arises in solutions 0{t) of rigid-body dynamics in the 
Lagrangian representation, where Ii denotes one princi- 
pal component of the inertia tensor for a symmetric top 
(/i = h 1^ h) and the effective potential V(cos^) con- 
tains terms associated with conserved angular momenta 
associated with the ignorable Eulerian angles and (f. 

Exact analytical solutions for these integrals exist 
only for certain potentials, in which cases the inversions 
t{x) x{t), 0{s) s{0), and t{0) 0{t) can be ex- 
pressed in terms of known functions. For example, exact 
solutions of the time integral ([T]) exist in terms of trigono- 
metric (or singly-periodic) functions when the potential 
U{x) is a quadratic polynomial in x. Trigonometric so- 
lutions of the orbit integral ([2]), on the other hand, exist 
for the Kepler problem U{r) = —k/r and the isotropic 
harmonic oscillator U{r) = kr'^/2. The purpose of the 
present paper is to explore exact analytic solutions of the 
integrals ©-(O expressed in terms of doubly-periodic 
functions called elliptic functions.^ For example, exact 
solutions of the orbit integral ([2]) for the central poten- 
tial U{r) = k exist in terms of elliptic functions^ for 
n = ±6, ±4,1, and —3. The time integral ([3|), on the 
other hand, has a solution in terms of elliptic functions 
for the problem of the heavy symmetric top (of mass M) 
with one fixed point (located at a distance h from the 
center of mass), where V(cos^) = Mgh cos 6. 



A. Doubly-periodic elliptic functions 



arises in solutions r{6) = l/s{0) of central- force problems 
involving the motion of a (fictitious) particle of (reduced) 



A function F{z) is said to be doubly-periodic, with 
periods r] and r]^ (where the complex- valued ratio r]'/r] 
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FIG. 1: Unit cell for a doubly-periodic elliptic function with 
periods 77 and ?7^ A doubly-periodic function F{z) satisfies the 
property F( 2; + m?7 + n?7^) = ^(z) for m, n = 0, ±1, ±2, ... 
(but not m = = n). 



has a positive-definite imaginary part), if 

F{z^mr]^nr]') = F{z), 

for m^n = 0, ±1,±2,... (but not m = = n). Figure 
[1] shows the unit cell (or fundamental period parallelo- 
gram) with corners at 2: = 0, 77, and r] -\- rj' in the 
complex plane. We note that, in the limit \r]'\ 00, 
the function F{z) becomes singly-periodic with period r]. 
Elliptic functions are doubly-periodic functions with 2 
simple zeroes per unit cell and either a second-order pole 
(Weierstrass elliptic function) or two first-order poles ( Ja- 
cobi elliptic function). Note that there are no elliptic 
functions of first order and that there are no multiply- 
periodic functions with more than two periods.^ 

Elliptic functions ^(x; a) are defined as solutions of the 
nonlinear ordinary differential equation 



— = a^y ^ asy ^ a2y ^ aiy ^ ao, 



where a = (ao, ai, a4) are constant coefficients. This 
equation can be formally solved by finding the inverse 
function 



x{y;a) = xo± 



f 



ds 



2/0 (a) Vci4 + a3 + a2 52 



• ai 5 + ao 



where ^0(^1) is a root of the quartic polynomial a4y^ + 
y^ 0.2 y'^ cii y -\- ao and x{yo; a) = xq. Jacobi elliptic 
functions are defined in terms of the quartic polynomial 
(1 — ?/^) (1 — m?/^), where m is a positive constant, while 
Weierstrass elliptic functions are defined in terms of the 
cubic polynomial 4y^ — g2y — gs^ where ^2 and gs are 
constants. There is a connection between the Jacobi and 
Weierstrass elliptic functions^ that will be exploited later 
in Sec. IIIIBI (see also Appendices [A] andlB]). 



B. Organization 

The remainder of the paper is organized as follows. 
In Sec. ini we present the Jacobi elliptic functions and 
discuss the Seiffert spherical spiral^^^ as a mathematical 
introduction of their doubly-periodic nature. Next, we 
discuss exact solutions to the physical problems of (i) 
the periodic motion in a quartic potential, (ii) the pla- 
nar pendulum, and (iii) Euler's equations for a force-free 
asymmetric top. In Sec. IIIH we present the Weierstrass 
elliptic functions and discuss exact solutions to the phys- 
ical problems of (i) the motion in a cubic potential, (ii) 
the planar pendulum (which demonstrates the connec- 
tion between the Jacobi and Weierstrass elliptic func- 
tions), (iii) the spherical pendulum, and (iv) the motion 
of a heavy symmetric top with one fixed point. In Sec. lIVi 
we discuss one interesting application of elliptic functions 
in terms of the travelling- wave solutions of nonlinear par- 
tial differential equations. We summarize our work in 
Sec. |V] and present mathematical details in Appendices 
[Al and [B1 Lastly, we note that the notation used for 
the Jacobi and Weierstrass elliptic functions in our pa- 
per is partly based on the standard material presented 
elsewhere jii^ 



II. JACOBI ELLIPTIC FUNCTIONS 

We begin our introduction of elliptic functions with the 
more familiar Jacobi elliptic functions. The Jacobi ellip- 
tic function sn{z \ m) is defined in terms of the inverse- 
function formula 
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sin Lp 



m sm t 
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- y'^) (1 -my^) 

= sn~^ {sin (f \ m) ^ (4) 

where the modulus m is a positive number and the am- 
plitude (p varies from to 27r. From this definition, we 
easily check that sn~^(sin(/:? | 0) = sin~^(sin (/:?) = cp. The 
solution to the differential equation 



(5) 



is expressed in terms of the Jacobi elliptic function 

sn{z\m) (for m < 1), 



m sn [rn^l^ z \ 



m 



(for m > 1). 



(6) 

By using the transformation y = sin the Jacobi differ- 
ential equation ([5]) becomes 



— — = 1 — m sm (z?, 
dz 



(7) 



FIG. 2: Plots of (a) sn{z\m) and (b) —isn{iz\m) form = 1/16 
showing the real and imaginary periods 4 K{m) and 4iK' (m) . 



FIG. 4: Plots of sn(2;|m), cn(2;|m), and dn(2;|m) from z = 
to 4:K(m) for m = (0.9)^ 



K' = K(m') 



^K= Kim) 



FIG. 3: Plots of the quarter periods K 
K{m') =K{l-m). 



K(m) and 



sm 



[sn(z|m)] 



and the solution to this equation is (p{z) 
for m < 1. 

The function sn{z\m) has a purely-real period 4i^, 
where the quarter-period K is defined as 



K 



K{m) 



.7r/2 

Jo 



dO 
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m sm 



and a purely-imaginary period 4:iK', where the quarter- 
period K' is defined as (with the complementary modulus 



1 — m) 



p7T/2 

IK' = iK{m') = i / 
Jo 



dO 
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m' sm 



Figure [2] shows plots of snz and —isn{iz) for m = 1/16, 
which exhibit both a real period and an imaginary pe- 
riod. Note that, while the Jacobi elliptic function sn2: 
alternates between —1 and +1 for real values of z (with 
zeroes at 2nK)^ it also exhibits singularities for imag- 
inary values of z at (2n + l)iK' (n = 0,1,...). Fur- 
thermore, as m ^ (and m' 1), we find K ^ 7:/2 
(or AK ^ 27t) and \K'\ oo (see Figure [3j), and so 
snz ^ sin z becomes singly-periodic. 



The additional Jacobi elliptic functions cn{z \ m) and 
dn(z I m) are defined as 



J cn 



dy 



{z\m) y/{l-y^){m' ^my^) 
dy 



(10) 



dniz\m) V(^-y'^)iy^ -^0 
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with the properties cnz 
dn{z\m) 



cn(2:|m) 



cos (/p, dn z 

\/\ — m sin^ and sn^z + cn^2: = 1 



dn z + msn'^^;. The Jacobi elliptic functions cnz and 
dn z are also doubly-periodic with periods 4 K and 4z K' 
(see Figure H]). 

The following properties are useful. First, we find the 
limits: 



and 



sn(z|0) 
cn(z|0) 
dn(z|0) 



/ sn(z|l) 
cn(z|l) 
Vdn(z|l) 




(12) 



tanh z 
sechz 
sech z 



(13) 



Next, we find the derivatives with respect to the argu- 
ment Z'. 



sn'(2:|m) = cn(z|m) dn(z|m) 
cn'(z|m) = — sn(2:|m) dn(z|m) 
dn'(z|m) = — m cn{z\m) sn{z\m) 

and, if m > 1, the identities: 

sn(z|m) = m~^^'^ sn{m^^'^ z\m~^) 
cn{z\m) = dn(m"'^/^ ^Im""*-) 
dn(2;|m) = cn{m^^'^ z\m~^) 



(14) 



(15) 



We now turn our attention to solving mathematical and 
physical problems using the Jacobi elliptic functions (|4|) 
and CT-(pT]). 



FIG. 5: Seiffert Spherical Spiral: Plot of the vector {[T7|) on 
the surface of the unit sphere for k = 0.15 from s = to (a) 
s = 2K, (b) s = 4K, (c) s = 6K, and (d) s = 8K. 



FIG. 6: Two different views of the Seiffert Spherical Spiral 
for k - 0.95 from s = to s = AK. 



A. Seiffert spherical spiral 

A simple example that clearly demonstrates the pe- 
riodicity of the Jacobi elliptic functions snz and cnz 
is given by the Seiffert spherical spiral^^'^ defined as a 
periodic curve on the unit sphere and constructed as 
follows. First, we use the cylindrical metric ds'^ — 



dz' 



vith 



yj\ — p^ and the az- 



imuthal angle Lp(^s) = ks is parametrized by the arc 
length s (assuming that the initial point of the curve 
is p = 0, (/p = 0, and z = 1). Hence, we readily find 
ds'^ = [(1 - (1 - k'^ p'^)]-^dp'^, which leads to (with 
< m = A:^ < 1) 



= r ^ 



dy 



= sn ^(p|m), 



V^l-t/2)(l-mt/2) 
and, thus, we obtain the Jacobi elliptic solutions 
p{s) = sn(5|m) 



z{s) = \/l — p'^{s) = cn(5|m) 



(16) 



The case m > 1 is handled with the identities (p!5|) . The 
Seiffert spherical spiral is generated by plotting on the 
unit sphere the path of the unit vector 

7{s) = sn(s|/c^) [ cos(/c5) X + sin(/c5)y] + cn{s\k'^)z 

(17) 

as a function of s (see Figure [5]). Note that at each value 
AnK {n = 1,2, ...), the orbit returns to the initial point 
at p = and z = 1. Figure [6] shows the complex periodic 
nature of the Seiffert-spiral orbit for k = 0.95. 



B. Motion in a quart ic potential 

Our next example looks at particle orbits in the (di- 
mensionless) quartic potential U{x) = 1 — + x^/16 




FIG. 7: Quartic potential U{x) = 1 — + showing 
orbits with ^ > 1, ^ = 1 (separatrix) and E < 1. 



shown in Figure [3 Here, the turning points for E = = 

U{x) are 



±2 vTTe (for e > 1) ' 
and ± 78 (for e = 1) 



> . 



(18) 



±2 vTTe (for e < 1) ^ 

Each orbit is solved in terms of the integral ([1]) using the 
initial condition = 2 \/l + e with the initial velocity 
xq < 0: 



t{x) 



L 



dy 



= - I 

\/e Jo 



2Vi+e V2(e2-1) + y2(l_y2/8) 
r VSdy 
Ayi+i v'[4(e+l)-y2][y2 + 4(e-i 



(19) 



y/l — m sin^ cp 

where m = (1 + e)/2e while we used the trigonometric 
substitution y = 2 yi + e cos (p with 



$(x) = cos ^ 



_2vTT^_ 



(20) 
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FIG. 8: Phase portait for orbits (|2T]) - (|23)) of the quartic po- 
tential U{x) = 1 - + x^/16 for (a) e > 1, (b) e = 1 
(separatrix) , and (c) e < 1. 



to obtain the last expression in Eq. (p!9|) . The Jacobi 
elliptic solutions obtained from Eq. ([19]) are shown in 
Figure [8] for the orbit (a), with e > 1, the separatrix 
orbit (b), with e = 1, and the orbit (c), with e < 1. 
For e > 1 (i.e., m < 1), we use Eq. ^ to find 



sin$(x) = sn(v^t| 



V 4(l + e)' 
which yields the phase-portrait coordinates {x^x): 



v{t) = 2^/TTe cn(v/et|m) 



x{t) = - 2 Y^e(l + e) sn{y/e t\m) dn{y/e t\m) 

'(21) 

where the velocity x{t) is obtained by using Eq. ([Hj). 
For e = 1 (i.e., the separatrix orbit with m = 1), the 
phase-portrait coordinates become 



x{t) = \/8 secht 

x{t) = — secht tanht 



(22) 



where the limits (fT3|) were applied to Eq. ([2T|) . Lastly, 
for e < 1 (i.e., m > 1), we apply the relations (p!5|) on 
Eq. ([2lD to obtain 



(23) 



x{t) = 2vTTe dn(r|m-i) 

x{t) = — >/8 e sn(r |m~-^) cn(r |m~-^) 



where r = t v'(l + e)/2. The orbits (f2T]) - ([23l) are com- 
bined to yield the phase portrait for the quartic potential 
shown in Figure [51 



C. Planar pendulum 

As a first physical example, we consider the well-known 
planar pendulum with a normalized energy ccOq {cOq = 



FIG. 9: Phase portait for pendulum orbits: (a) bounded (os- 
cillation) orbits (m < 1) and (b) unbounded (libration) orbits 
(m > 1). The separatrix orbit (dashed line) separates oscilla- 
tions from librations (which are periodic since — tt is identical 
to tt). 



g/L) expressed as 



1 • 



r + cjg(l -cos(9) 
2lo^ (^'2 + sinV) , 



(24) 



where (f = 0/2 and (f'{r) = coQ^cp. The normalized time 
(assuming that (/p = at r = 0) is thus expressed as 



pep 

^0 ^ 



(e/2) -sin^ 



(25) 



From Eq. ([24|) . the dimensionless energy e can either be 
(I) < e < 2 (i.e., (p' vanishes along the orbit) or (II) 
e > 2 (i.e., (p^ does not vanish along the orbit). In case 
(I), we set e = 2 sin^ a and smcf) = sin a sinx in Eq. ([25]) 
to obtain^ 



arcsmi m 



L 



dx 



\/l — m sin^ X 



where the modulus is m = sin^ a = e/2. Hence, the 
solution for case (I) is expressed in terms of the Jacobi 
elliptic function 



sm 



m"*"/^ sn(r | m) 



(26) 



In the limit m <C 1 (i.e., e <C 2), we use the limit ([T2]) to 
obtain the simple harmonic solution 9{t) = \/2e sinr. 
In case (II), we have = 2/e < 1 and obtain 



[2 r 

V 6 io TT 



m sm (p 
m~"^/^ sn~"^(sin | m~^). 
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Hence, the solution for case (II) is expressed in terms of 
the Jacobi ehiptic function 



sn(m"^/^r I m ^) 



(27) 



which fohows from property (p!5|) . In the hmit m ^ 1, 
both solutions ([26|) and ([27j) coincide with the separatrix 
solution 



tanhr, 



(28) 



which is expressed in terms of singly-periodic (hyper- 
bolic) trigonometric functions (with imaginary period). 
Since cp ±7r/2 (^ ^ ±7r) as r ^ ±oo, the period 
of the pendulum on the separatrix orbit is infinite. The 
pendulum orbits ([26|) - ([28|) are shown in Figure [9l 

We will return to this example in Sec. IIIIBI where we 
solve the problem of the planar pendulum in terms of the 
Weierstrass elliptic function and establish a connection 
between the Jacobi and Weierstrass elliptic functions. 



D. Force- free asymmetric top 

As a second physical example, we consider the Euler 
equations for a force-free asymmetric top (with principal 
moments of inertia Ii > I2 > ^3)*^ 



hibi = {I2 - I3) ^2^3 
I2UJ2 -(71-/3)^1^3 

^3^3 = {h - h) 0J\ijJ2 



(29) 



where the angular velocity = c^i 1 + cc;2 2 + cc;3 3 is de- 
composed in terms of its components along the principal 
axes of inertia. The conservation laws of kinetic energy 



1 



1 



I A (30) 



and (squared) angular momentum 
e = IIloI + Ilul + Ilul 



t2 q2 

1q \Lq^ 



(31) 



are expressed in terms of the parameters Iq = £^/(2/t) 
and Qq =2k/£. These conservation laws can be used to 
introduce the following definitions 



wi(t) 



W2(r) 



'^3(t) 



V h {h - h) 



lo (Jo - h) 

V h ih - h) 

Wo) y{T), 

l loih-Io) 

V ^3 (/l - h) 



^0 y(r) 



(32) 



(33) 



(34) 




FIG. 10: Plots of {001,002,003) at different times: (a) r 
(b) r = 2K, (c) r = 3K, and (d) r = 4K. 



where r = [(/i — /a) Vti ^2)] t is the dimensionless 

time and the modulus m is defined as 



m(/o) 



(/Q-/3) (/1-/2) 
(/2-/3) (/1-/0)' 



(35) 



By requiring that the modulus m be positive, the pa- 
rameter Iq introduced in Eqs. ([3Q|) - ([3T]) must satisfy 
Is < Iq < Ii and, hence, < m(/o) < 1 for I3 < Iq < I2 
and m{Io) > 1 for I2 < Iq < h (with m ^ 00 as 
/o^/i). 

When we substitute these expressions in the Eu- 
ler equation ([29|) for 002 = 1^2 ^(t), we easily ob- 
tain the dimensionless Jacobi differential equation ([5]), 
which can now be integrated, with the initial conditions 
(c^i(0),cj2(0),c^3(0)) = (-1^1,0,(^3), to yield^ 

(cji, 002^ ^3) = (— l^icnr, 1^2 snr, 1^3 dnr). (36) 

These solutions are shown in Figure [10] where the AK- 
periodicity is clearly observed. The solution of this prob- 
lem is thus very elegantly expressed in terms of the Jacobi 
elliptic functions (sn, cn, dn). 

The separatrix solution (m = 1) corresponds to the 
case when Iq = I2 for which 0.2 = ^0, so that the sepa- 
ratrix solution is 



cji(r) 
^3(t) 



hjh-h) 
V h - h) 
$7o tanh t, 



1^0 sechr, 



I hih-h) 
h {h - h) 



sechr. 



Lastly, a symmetric top (/i = I2 7^ ^3) corresponds 
to the limit m = (and Q2 = ^i)- The Jacobian 
solution ([36]) for a symmetric top therefore becomes 
(cji, 002, ^3) = (— ^1 COST, Oi sinr, 1^3), where we 
have used the identities ([12]) and r = (1 — /3//1) 1^3 t is 
the dimensionless time. 

The examples presented in this Section dealt with peri- 
odic orbital motions described in terms of Jacobi elliptic 



7 




= 4(5-ei)(5-e2) (5-63). (37) 



Here, (ei, 62, 63) denote the roots of the cubic polynomial 
45^ — g2S — gs (such that ei + 62 + 63 = 0), where the 
invariants ^2 and gs are defined in terms of the cubic 
roots as^^^ 



A<0 \ A>0 \ A<0 



FIG. 11: Cubic roots (61,62,63) as a function of e = 
(3/5'2)^^^^3 with fixed value ^2, where a = y^^g2jl2 and 
^ — 92 ~ • The three roots satisfy 61+62 + 63 = 0. 
The regions (a)-(c) along the e-axis correspond to segments 
of the path shown in Fig. 1121 



92 



93 



-4 (ei 62 + 62 63 + 63 61) ' 

2 {el + 6^ + el) 



4 61 62 63 



(38) 



functions (sn, cn, dn) with arguments evaluated along the 
real axis (i.e., the two simple-pole singularities of the Ja- 
cobi elliptic functions are located on the imaginary axis 
of the Jacobi fundamental period-parallelogram). In the 
next Section, we consider the Weierstrass elliptic func- 
tion, which has a double-pole singularity on the real axis 
of the fundamental period-parallelogram. 

III. WEIERSTRASS ELLIPTIC FUNCTIONS 

The Weierstrass elliptic function p{z]g2^gs) is defined 
as the solution of the differential equation 

{ds/dzf = 4s^ - g2S - gs 



and A = ~ ^7 ^3 is the modular discriminant. Since 
physical values for the constants ^2 and gs are always 
real (and g2 > 0), then either all three roots are real or 
one root (say 6^) is real and we have a conjugate pair 
of complex roots (65, 6^) with Re(65) = —6^/2. The ap- 
plications of Weierstrass elliptic functions are analyzed 
in terms of four different cases based on the signs of 
(^3, A) = [(-,-),(-,+),(+,-),(+,+)], with two spe- 
cial cases (^3 7^ 0, A = 0) and (^3 = 0, A > 0). 

In general, the roots (61, 62, 63) of the cubic polynomial 
on the right side of Eq. (|37|) can be expressed in terms of 
the parameters a = ^/g2/l2 and e = (8/^2)^^^ gs as 









1 = a 


yes) 





(e + Ve^ - 1) 



-1/3 



Itt/S 
t/3 



1/3 



oi7r/3 



(39) 



and the discriminant is 

A = gi- 27 gl = gl{l-e'). (40) 

These cubic roots are shown in Figure [TTl as a function of 
e for a = 1/2 (i.e., g2 = 3); the polynomial 45^—^25 — ^3 
is positive (and ds/dz is real) to the left of the curve 
and negative (and ds/dz is imaginary) to the right of the 
curve. The roots can be further parametrized by writing 
e = — cos with e + — 1 = — exp(— i (p) and 



SO that the three roots (|39|) become 

61 \ / cos[{(p - 7r)/3] \ 

62 = 2a cos[((/:) + 7r)/3] . (41) 

63 y y - cos((/p/3) J 

The three cubic roots are connected smoothly as shown 
in Figure [11] by following the path in the complex (/p-plane 
shown in Fig. [121 Here, for e < — 1, the imaginary phase 
(p = —ill; (with > 0) yields 61 = a — z6 = 63 (with 
b > 0) and 63 = —2a< —1; for —1 < e < 1, the real 
phase = (j) (with < (j) < tt) yields 61 > 62 > 63; and 
for e > 1, the complex phase cp = tt -\- iip (with > 0) 
yields 61 = 2a > 1 and 62 = — a — i6 = 63 (with b > 0). 
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(a) 



(b) 



(c) 



Re(cp) 



(a) 


(b) 


Im(ro') 






V ^ Im(Q') 


0)-^ 





FIG. 12: Weierstrass path in complex Lp plane: (a) e < — 1, 
(b) -1 < e < 1, and (c) e > 1. 





FIG. 13: Plots of (a) > and (b) p(iz) < for ^2 = 3 

and — 0.5 (with e = 0.5) showing (a) the real period 2 a; 
and (b) the imaginary period defined, respectively, by 
Eqs. (14^ and (l43l). 



Note that 62 = (and e\ = >/3a = — 63) for ^3 = (i.e., 
(\) = 7r/2); this case is called the lemniscatic case.^ 

Figure [13] shows that, for < e < 1 (i.e., A > 0), 
p{z) has different periods 2uj and 2uj^ along the real and 
imaginary axes, respectively, with the half-periods co and 
uj^ defined as 



r 



ds 



(42) 



y^4s^ - 92 8- g^' 

^'{92.9^) = ^ r -7fTT—=^' 
J -00 ^/\4s^ - 92S - 93] 

The plots of uj{g2, 93) and cc;'(^2, ^3) are shown in Fig. [TH 
for ^2 = 3 as functions of 93. Note that for ^3 = 
(with ei = — 63 and 62 = 0), we find that uo' = iuo while 
\ui)'\ approaches infinity as ^3 approaches one. Explicit 
calculations of the half periods {001^002^003) in terms of 00 
and oj' are given in Appendix B. 

For e > 1 (i.e., A < 0), on the other hand, p{z) has 
different periods 2 Vt and 2 Vt' along the real and imagi- 
nary axes, respectively, with the half-periods and Q.' 



FIG. 14: Plots of (a) cj and lui{uj') for < ^3 < 1 and ^2=3 
(i.e., A > 0) and (b) Q and Im(r2^) for ^3 > 1 and ^'2 = 3 (i.e., 
A < 0). Note that u;'(^2,0) = iu;(^2,0), Q(^2, 1) = oo{g2, 1), 
and both (Q, Q') decrease to zero as ^3 becomes infinite. 



defined as 



^{92,93) 



r 



ds 



(44) 



4:S^ - 92 s - 93' 

f^^ ds 
'{92.93) = I / . 3 (45) 
J -00 \/\'^s^ - 92S - 93\ 



The plots of 1](^2, ^3) and Vt'{g2^g3) are shown in Fig. [Ml 
for ^2 = 3 as functions of ^3. Note that oo{g2.1) = 
Vt{g2^1)^ \^'\ approaches infinity as ^3 approaches one, 
and that both and Vt' approach zero as ^3 approaches 
infinity. 

Figure [15] shows the integration path used to calculate 
the half-periods {00, 00') for A > and {Vt, Vt') for A < 
0. When the integration path includes a segment where 
the polynomial At^ —92^ — 93 is positive (solid line in 
Fig. [T5|) , the contribution from this segment is a real- 
valued number. On the other hand, when the integration 
path includes a segment where the polynomial At^ —g2 t — 
g3 is negative (dashed line in Fig. [T5|), the contribution 
from this segment is an imaginary number. 

Table [I] shows the cubic roots = (61,62,63), defined 
by Eq. (|^T]), and the half periods ooi = {001^002.003)^ de- 
fined as^ 



i{92.93) / 



ds 



^ ^45^ - g2S- g3 
ds 



i 2 ^/{s- ei){s - e2){s - 63) 



(46) 



The cubic roots and half periods satisfy the following 
properties: 



P{^i) 

p{z^uji) = ei + 
p{z^2uji) = p{z), 



{Cj - Cj) {Cj - 6fc) 
P{z) - 6i 



(47) 
(48) 
(49) 



where i 7^ j 7^ A: so that p{oOi -\- ooj) = Ck- 



9 




Figure [16] shows the plots of p{z -\- UJ2) and p{z -\- ujs) 
for one complete period from 2: = to 2a;i, which clearly 
satisfies the identities (|47|) and (|49|) . The singular behav- 
ior exhibited by Eq. (|48|) at z = =b cji, on the other hand, 
is shown in Fig. [13] for i = 1 and j^k ^ 1. 



FIG. 15: Integration paths used to calculate the half-periods 
{uj, uj') for A > and (Q, Vl') for A < 0. Solid (dashed) 
lines indicate where the polynomial 4t^ — g2t — is positive 
(negative). 




-1 



FIG. 16: Plots of (a) p{z-\-uj2) and (b) ^(z+cjs) for ^2 = 3 and 
^3 = 0.5 (with e = 0.5) over one complete period from to 
2a;i. Note that p{<jOj) — Cj for j = 2 or 3 and p{ijJi-\-ix}j) — efc, 
for i = 1 and (j, k) = (2, 3) or (3, 2). 



TABLE I: Cubic roots (61,62,63) and half periods {001^002^002,) for 5^2 = 3 (a = 1/2). 



(^3 , A) 


(-, -) 


(-,0) 


(-, +) 


(0, +) 


(+, +) 


(+,0) 


(+, -) 


61 


a — ib 


1/2 


d>0 


^/2 


c> 


2 


2a > 1 


62 


a-\- ib 


1/2 


c-d>0 





d-c<0 


-1 


— a — ib 


63 


-2a< -1 


-1 


-c < 


-VS/2 


-d<0 


-1 


— a-\- ib 


001 


+2Q/2 


00 


kl 


3-^^^K{l/2) 


00 


Tt/a/S 


Q 


002 


- +iQ/2 


— 00 


icj — \oj^\ 


[3-^/^i^(l/2)](-l+z) 


— 00 — 00' 


— i 00 


-Q/2 - Q' 


OOz 




—iir/V^ 


—ioo 


-i3-^/^i^^(l/2) 


00' 


i 00 


-Q/2 + Q' 



The Weierstrass elliptic function p{z; ^2,^3) obeys the 
homogeneity relation-*^ 

p{Xz]X~^ g2,X~^ gs) = A"^ p(2:;^2,^3), (50) 

where A 7^ 0. By choosing A = —1, for example, we 
readily verify that the Weierstrass elliptic function has 



even parity, i.e., p{-z;g2,g3) = p{z;g2,g3)' On the other 
hand, for A = i, we find that the half-period assignments 
for ^3 < in Table [T] are based on the relation 

p(^;^2,^3) = - p{iz;g2, IgsD- (51) 
For example, for —1 < ^3 < (and A > 0), we find for 
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TABLE II: Bounded and unbounded orbits in a cubic potential (see Figs. [T7land ll8|) . 



Orbit 


Energy 


Time Range 


Constant 7 


Period 


Turning Point (s) 


(a) 


E > 2/3 


-iQ < t < iQ 




Unbounded 


63 < 


(b) 


0<E <2/3 


0<t< 2\uj'\ 


— ioo 


2\co'\ 


63 < 62 


(c) 


0<E <2/3 


-\uj'\ <t<\uj'\ 


ki 


Unbounded 


61 


(d) 


-2/3 < ^ < 


<t <2oo 


a.' 


2oo 


63 < 62 


(e) 


-2/3 < ^ < 


— UJ < t < (jJ 




Unbounded 


61 


(f) 


E < -2/3 


-Q < t < Q 




Unbounded 


61 




(c) 





xX 


(e) 






X,(ft — 


(t) 



convention^ for the case (^3, A) = (+,—), the root 62 
is real (while = 63) and the corresponding half-period 
002 is also real (in contrast to the convention adopted in 
Table H]). The connections (|52|) shown in Table [J are sim- 
ply based on the smooth dependence of the cubic roots on 
the single parameter e (for fixed ^2). These connections 
will enable us to describe consistent orbital dynamics in 
several problems in classical mechanics. 



A. Motion in a cubic potential 



FIG. 17: Cubic-potential energy levels E — x — x^/3 showing 
orbits (a) E > 2/3 (unbounded orbits; e < —1), (b) and (c) 
< £^ < 2/3 (bounded and unbounded orbits; — 1 < e < 0), 
(d) and (e) — 2/3 < ^ < (bounded and unbounded orbits; 
< e < 1), and (f) < -2/3 (unbounded orbit; e > 1). 



p(^i;^2,^3): 

p(l^'l;^2,- |^3|) 



p(^';^2, |^3|) 



i-d) 



which corresponds exactly to ei = d found in Table U for 
the case (^3, A) = (-,+). 

In general, the connections between the half-periods 

{iJi^iJ2i^t) ~^ {^1 1^2 ^^3) cubic roots 

(e^, 62", 63") (ej", 62 J 63 ) as gs changes sign from pos- 
itive (+) to negative (— ) are found in Table [J to be 



(^1,^2 ^^3) = (-^^3".-^^2".-^^^) 

(61,62,63) = (— 63", — 62", — e^) 



> . 



(52) 



Once again, these connections follow a non-standard 
convention. For example, according to the standard 



As a first example of problems solved by the Weier- 
strass elliptic function p, we consider particle orbits in a 
(dimensionless) cubic potential U{x) = x — Here, 
the cubic-potential orbits x{t) are solutions of the differ- 
ential equation 



x' 



2 E 



- {x — xi) {x — X2) {x ■ 

o 



X3), (53) 

and the turning points {xi^X2^xs) are shown in Figure 
[TTl (with xi -\- X2 -\- xs = 0). By writing x{t) = 6s{t), 
Eq. ([53|) is transformed into the standard Weierstrass el- 
liptic equation ([37|) . where the invariants are g2 = 1/3 
and gs = -E/18, so that e = -3E/2. Note that 
bounded orbits exist only for — 1 < e < 1 (i.e., A > 0). 

The cubic-potential solution is therefore x{t) = 6 p{t-\- 
7), where the constant 7 is determined from the initial 
condition x{0). In Fig. [181 the orbits (a)-(f) are shown 
with initial conditions identified by a circle and a quali- 
tative description of these orbits is summarized in Table 
im Note that the turning points Xi = Ge^ {i = 1,2,3) 
are simply related to the standard cubic roots e^. Lastly, 
the separatrix solution is obtained from orbit (b) as E 
approaches 2/3 and the period 2 becomes infinite. 



Note that the imaginary time range for orbit (a) takes as the single (energy) parameter e is varied, 
into account the relation ([5T]) since ^3 < for this orbit. 
In addition, the connections ([52|) allow us to describe the 
orbits (a)-(f) in Figs. [l71 and [18] (and Table[IIl) smoothly 
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FIG. 18: Plots of x(t) versus x(t) for cubic potential ([53|) 
shows bounded and unbounded orbits: Orbit (a) E > 2/3 
(e < -1 and A < 0); orbits (b)-(e) -2/3 < ^ < 2/3 (-1 < 
e < 1 and A > 0); and orbit (f) E < -2/3 (e > 1 and A < 0). 
The dotted lines are the bounded and unbounded separatrix 
orbits for E = 2/3 and circles denote particle positions at 
t = 0. 



B. Planar pendulum 

We now return to the planar pendulum problem of 
Sec. Ill CI where we write z = 1 — cos 6 (i.e., < z < 2) 
and transform Eq. into the cubic-potential equation 

{z'f = 2z{2-z){e-z), (54) 

with roots at z = 0, 2 and e. When e < 2, the motion 
is periodic between z = and z = e, while the motion 
is periodic between z = and z = 2 for e > 2. We re- 
cover the standard Weierstrass differential equation ([37|) 
by setting 

z{r) = 2p(r + 7) + /i, (55) 

where /i = (e + 2)/3 and the constant 7 is determined 
from the initial condition z{0). 

The root corresponding to 2: = is labeled Cc = — 
the root corresponding to z = 2 is labeled 65 = 1 — /i/2, 
and the root corresponding to z = e is labeled = /i — 1 
and we easily verify that + 65 + Cc = (see Fig. [T9|) . 
The Weierstrass invariants ^2 and gs are 

92 = 1 + 3 (/i - 1)^ and gs = /i (/i - 1) (/i - 2), 

and the modular discriminant is A = (2 — e)^ > 0. 

The planar pendulum is discussed in terms of 4 cases 
labeled (a)-(d) in Fig. [191 For cases (a) and (b), where 
2/3 < /i < 4/3 (i.e., < e < 2), we find 63 = - /i/2 < 
e2=/i — l<ei = l — /i/2, so that = {ei — e^)^!^ = 1 
and m = {62- es)/{ei - 63) = (3/i - 2)/2 = e/2 < 1. 
For cases (c) and (d), where /i > 4/3 (i.e., e > 2), we 
find 63 = — /i/2 < 62 = 1 — /i/2 < ei = /i — 1, so that 




(a) Cb) (c) (d) 



12 4 8 



FIG. 19: Plots of the cubic roots (ea,eb,ec) as functions of 
/i = (e+2)/3. The cases (a)-(d) are discussed in the text. Note 
that for cases (a) and (b), or e < 2, we find Cc < Ca < eb, 
while for cases (c) and (d), or e > 2 we find Cc < et < ea- The 
bounded motion of the planar pendulum ( — 1 < ^ < 1) occurs 
between the two lowest cubic roots: Cc < Ca (for e < 2) or 
Cc < 66 (for e > 2). 




FIG. 20: Plot of the Weierstrass invariant ^3 as a function 
of e. For case (a), ^3 > and A > 0; for cases (b) and (c), 
^3 < and A > 0; and for case (d) ^3 > and A > 0. 

= (ei - 63)'/' = (e/2)'/' and m = (62 - es)/{ei - 63) = 
2/e < 1. Figure [2Ql shows a plot of g^ as a function of 
the parameter e, which can be used with the information 
presented in Table [J to describe the motion of the planar 
pendulum in terms of the Weierstrass elliptic function. 

We first consider case (a), where < e < 1 (i.e., 
2/3 < /i < 1 and ^3 > 0), the periodic motion is bounded 
between 63 = —/i/2 (i.e., z = 0) and 62 = /i — 1 < (i.e., 
z = e). Using the initial condition 2;(0) = 0, we find that 
pin) — = 63 which implies that 7 = 0;' (see uo^ in 

Table [J for ^3 > and A > 0). The Weierstrass solution 
of the planar pendulum for < e < 1 is therefore 

z{t) = 2p{T^u') + /i, (56) 

with the period of oscillation 2uj. For the case (b), where 
1 < e < 2 (i.e., 1 < /i < 4/3 and g^ < 0), the periodic 
motion is bounded between 63 = —/i/2 (i.e., z = 0) and 
62 = /i — 1 < (i.e., z = e). Using the initial condition 
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z{0) — 0, we find that ^(7) = — /i/2 = 63 which implies 
that J = — icj (see uos in Table [J for ^3 < and A > 0). 
The Weierstrass solution of the planar pendulum for 1 < 
e < 2 is 



z{t) = 2 p(r - iuj) + /i, 



(57) 



with the period of oscillation 2 \u;'\. As expected, when 
e ^ 2 (i.e., A ^ and m ^ 1), the period 2 \u;'\ ap- 
proaches infinity as we approach the separatrix. For case 
(c), where 2 < e < 4 (i.e., 4/3 < /i < 2 and ^3 < 0), 
the periodic motion is bounded between 63 = —/i/2 (i.e., 
z = 0) and 62 = 1 — /i/2 (i.e., z = 2), with the period of 
oscillation 2|a;'|. Using the initial condition 2:(0) = 0, 
we find that ^(7) = —/i/2 = 63 which implies that 
J = —iu (see ujs in Table [J for ^3 < and A > 0) 
and thus the Weierstrass solution of the planar pendu- 
lum for 2 < e < 4 is again given by Eq. ([57j) . Note that 
the separatrix solution (e = 2) is represented by orbits 
(b) and (c) as \u;'\ 00. Lastly, for case (d), where e > 4 
(i.e., /i > 2 and ^3 > 0), the periodic motion is bounded 
between 63 = —/i/2 (i.e., z = 0) and 62 = 1 — /i/2 (i.e., 
z = 2). Using the initial condition z{0) = 0, we find that 
Pil) — = 63 which implies that 7 = cj' (see CJ3 in 

Table U for ^3 > and A > 0) and thus the Weierstrass 
solution of the planar pendulum for e > 4 is again given 
by Eq. ([561). 

We conclude our discussion of the planar pendulum by 
using the Jacobi and Weierstrass solution of this prob- 
lem to establish a general relation between these elliptic 
functions. First, we use the Jacobi elliptic solution ([26|) 
for the case e < 2 and find 



z{t) = 2 sin^ 



2 m sn^(/^ r | m). 



(58) 



where m = e/2 and hi = 1. By comparing Eqs. ([56|) and 
([581). we obtain the relation 



/i 
2 



; /i — 1 



sn 



1 



which is just one example of the general relation 

p(r + cj3) = 63 + (62 - 63) sn^(/^r I m). (59) 



We note that this relation validates the relation (|Bip be- 
tween the Weierstrassian half-period ui = u and the Ja- 
cobian quarter-period K = nuo. Next, using the Jacobi 
elliptic solution ([27|) for e > 2, we find 



(60) 



z{t) = 2 sn^ (^Ve/2r \ 2/e^ 



2 sn^ 



and thus we recover once again the relation ([59|) . 



C. Spherical pendulum 

Our second physical example is the spherical pendu- 
lum (whose pendulum bob moves on the surface of a unit 



sphere), which is described here in terms of the cylin- 
drical coordinates (p^cp^z)} The energy equation for a 
spherical pendulum of unit mass and unit length is 



1 



•2 I 2 -2 I ■2\ 



+ z, 



where e denotes the normalized total energy {u'^ = g). 
By substituting p{z) = \J\ — z^ (where — 1 < ^ < 1) and 
the angular momentum conservation law i = we 
obtain the differential equation 



i\2 



2{e-z){l-z'^) - 
2{z- zi) {z - Z2) {z - Z3), 



(61) 



Because 



where z\r) = z and 2:1 + ^2 + ^3 = 
the right side of Eq. ([6T]) is negative at z = ±1, the 
highest root zi (> 1 > Z2 > zs > — 1) of the cubic 
polynomial is greater than 1 and is, therefore, unphysical 
(since the cylindrical radius p then becomes imaginary). 
The periodic motion of the spherical pendulum is thus 
bounded between zs < z < Z2. 

The differential equation ([6T]) can be transformed into 
the standard differential equation ([37|) by setting z{r) = 
2 p{r + 7) + /i, where p = e/3, the constant 7 is deter- 
mined from the initial condition z{0), and the invariants 
g2 and gs are ^2 = 1 + 3 /i^ and gs = ^^/4 + p {p^ — 1). 
Note that, if z(0) = ^3 > — 1, then 7 = 0^3 (i.e., uo' 
when A > and Vt' when A < 0) and the solution of the 
spherical pendulum problem for the z-coordinate is 



z{t) = 2p(r + c^3) + /i = Vl - P^{r). (62) 

The motion is periodic with period 2a;i (i.e., uo when 
A > and ^ when A < 0) on the unit circle in the (p, z)- 
plane. At the half-period c^i, we find z{uoi) = 2p{uoi + 
^3) + /i = 2 p{(j02) + /i = 2^2 < 1 as expected. 

The solution for the azimuthal angle (/p(r) is obtained 
from the angular- momentum conservation law ^'{t) = 
ij P^(t\ which is integrated as-^ 



= '/'- 



ds 



[2p(5 + cc;3)+/i]2 
ds 



(63) 



4 Jo [p(<5 + ^3) - p(^)] \p{s + ujz) - p{\)\ ' 

where we used the initial condition (/p(0) = and the 
imaginary constants n and A are defined by the relations 
p(n) = — (1 + /i)/2 and p(A) = (1 — /i)/2 corresponding 
to z = — 1 < Z3 and z = + 1 > Z2 , respectively. These 
constants also yield the relations p' (n) = i£/2 = p'(A), 
obtained from the Weierstrass differential equation ([37|) 
for z = hz and A. These relations allow us to write Eq. ([63]) 
as 



- f 

2 7o 



ds 



P'(A) 



p(s + wa) - p{\) 
P'(^) 



p(s + W3) - p(k) J ' 



(64) 
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where we used the identity p(A) — p{hz) = 1. 

The integral can be solved exactly in terms of 
the quasi-periodic functions ({r) and cr(r) associated 
with the Weierstrass elliptic function p{r):^ p{r) = 
— d({T)/dr and ({r) = a' {t) / cf{t). Using the identity 



2C(A) 



ds 



In 



a{s-\y 
a(s + A), 



P'(A) 



p(s) - p(A) 



C(5-A) - C(5 + A) + 2C(A) 



we find the standard solution^ for the azimuthal motion 
of the spherical pendulum 



,2i(^(r) 



,2r [C(/^)-C(A)] 



cr(r + cja + A) cr(cc;3 - A) 
cr(cc;3 + A) cr(r + 6^3 - A) 

I 



cr(r + cc;3 - cr(cj3 + n) 



(65) 



where we can easily verify that (/?(0) = 0. This solu- 
tion can be studied numerically with Mathematica since 
its library of functions contains the Weierstrass elliptic 
function p(r) and its associated quasi-periodic functions 
C(r) and cr(r). 

We now simplify Eq. (|65|) by recognizing that, since 



J 



— 1 < ^3 ^ pij + ^3) < 2:2 < + 1 for real values of r, 
there must be imaginary numbers i a and i (5 (where a 
and (3 are real- valued constants) such that n = uo^, ^ i a 
and A = CJ3 + z /3. Using these substitutions, the solution 
(|65|) becomes 



,2r [C(/^)-C(A)] 



a(r + 2cj3 + i/?)cr(-i/^) 

Cr(2cJ3 + ^/?)cr(r-i/3) 



cr(r — z a) cr(2 cc;3 + i a) 
(t(— m) (T(r + 2 cc;3 + i a) 



(66) 



Next, using the identity^ cr(r + 200^) = — exp[2 7^3 (r + cc;3)] cr(r), where 7^3 = C(^3)7 we now find 



and 



(j(t ^2uo2.^ip) 



a(r-z/3) 



^-2i 773/3 ^( ^/^) = g-2i 773/3 



Cr(2cJ3 + ^/5) ^ Cr(i/5) 

where we used the fact that <j(r) is an odd function of r, so that i/3) = — /3). Hence, we obtain 



and 



(j(r + 2cj3 + i/3)a(-z/3) 
(7(2003 ^ i P) cr{T -ip) 

cr(r — ia) (j{2 6^3 + i a) 



-2 773r 



cr(— za) cr(r + 2 CJ3 + i a) 
We combine these relations to obtain the simplified solution 



(j{r — ia) 
_ a{r -\- ia) _ 



,2i(^(r) 



,2r [C(/^)-C(A)] 



cr(r + i cr(r — i a) 



_ cr(r — i /?) cr(r + z a) _ 

Lastly, we note that the ratio cr(r + i /3)/(j{r — i/3) must have unit modulus for real r- values and so we may write 



(67) 



In 



(T{r^i(3) 
a{r-ip) 



I / C,{t -\- i s) ds 

J-f3 

2i / Re[C(r + zs)] ds. 
Jo 
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Hence, the azimuthal angle finally becomes 



ip{r) = / Re[C{r ^ i s)] ds + ir[C{uJ3^iP) - C(^3 + ^Q^) 

J a 

= Re ij [C(r + i 5) + r p{u;s i s)] dsj , 



(68) 



r 



which is now expressed only in terms of the quasi-periodic 
function C(r) and C'('^) = ~p{'t)- After a full period 
2a;i, when the (p, z)-coordinates return to their initial 
values, the azimuthal angle has changed by an amount 
Acp = (p{t + 2a;i) — (p{r) expressed as 

A(p = 2uJi / p{uj3 + is)ds + 2rii{f3-a), (69) 

'J a 

where we used the identity^ C{^~^ s) = ({r s) -\- 

2r]i and r]i = C(^i). 



D. Heavy symmetric top with one fixed point 

Our third physical example is provided by the motion 
of a symmetric top (/i = h 7^ h) with one fixed point 
described in terms of the energy equation 



E 



Mgh cosO. 



[Pip -Pip cos i9)^ 



(70) 

where the total energy E of the symmetric top (with total 
mass M and principal moments of inertial h = I2 7^ 
/s), ws denotes the constant component of the angular 
velocity, and the angular momenta p^p and associated 
with the ignorable Eulerian angles (f and are constants 
of motion. By defining the dimensionless parameters e = 
{E - \hwl)/{Mgh), {a,b) = {p^/ hu, p^/ hu), where 
z/2 — Mgh/{2 Ii), the differential equation for u = cos^ 
is obtained from Eq. ([7Q|) as 

{u')^ = 4{l-u^){e-u)-{a-bu)^ 

= 4:{u - ui) {u - U2) {u - Us), (71) 

where a prime denotes a derivative with respect to the 
dimensionless time r = vt and < U2 < ui are the 
three roots of the cubic polynomial. Since the right side 
of Eq. ([TT]) is negative at = ±1, we conclude that 
— l<us<U2<l and Ui > 1 (which is unphysical for 
u = cos^). The physical motion is therefore periodic in 
and is bounded between = cos 6^3 and U2 = cos ^2 
(or ^2<^(r)<^3). 

By using the change of integration variable u = s -\- fi, 
where fi = (4e + 6^)/12, the differential equation ([TT]) 
becomes the standard differential equation (|37|) for the 
Weierstrass elliptic function, with 



g2 = 2 {2 - ab ^ 6/i^) , 
^3 = - 4 e + 2 /i (2 - a6) 



Hence, the solution is expressed in terms of the Weier- 
strass elliptic function as 



u{r) cos 6>(r) = p(r + 7) + /i. 



(72) 



where 7 is determined from the initial condition ^(0). 
Assuming that —l<us = es-\-fi<U2 = e2-\-/J^< 
1 < iii = ei + /i, we choose u{0) = es -\- fi = us (i.e., 
— 1 — /i < 63 < 1 — /i) so that J = co' and, hence, at 
the half-period r = cj, we find u{u;) = p{u; + co') -\- jj. = 
62 + /i = '^2- The solution for 6{r) is thus expressed as 
l9(r) = cos-^[p(r + cc;0 

Note that the sign of cp^ depends on the sign of a — 
b cos6>2 < a — b cos6 < a — b cos ^3. If a > 6 cos ^2 (or a < 
b cos ^3), cp' does not change sign as bounces between O2 
and Os and the motion in (f involves monotonic azimuthal 
precession. If a = 6 cos ^2 (or a = 6 cos 6^3), (p' vanishes at 
6 = 62 {oT 6 = Os) and the motion in (p exhibits a cusp at 
that angle (since both 0' and cp' vanish). If a < 6 cos ^2, 
(p' vanishes at an angle O2 < Oq < Os and the motion in 
ip exhibits retrograde motion between ^0 < < Os- 

The solution ([72|) for the Eulerian angle 0{t) can now 
be used to integrate the differential equations 



a — b cos 0{r) 

i 
2 



cos^ 0{r) 



(73) 



and, defining iIj' 



p(^) 
iws/i^ - b) - 



-^3) - p(A) 



X'(r) 



b — a cos 0{r) 
1 - cos2 (9(r) 



(74) 



P'(A) 



-^3) 



p(-^) 



• + CJ3) - p(A) J ' 



for the remaining Euler angles, where p{hz) = 1 — /i 
and p{X) = — (1 + /i), with p'{tv) = i {a — b) and 
p'{X) = i {a -\- b). Note that since Eq. ([73|) is the same 
as Eq. (|64j) . its solution is identical to Eq. (|68|) (even if 
the constants Hi and A are different). This same solution 
can also be applied to the solution for Eq. ([74|) . where 
Eq. ([73|) is transformed into Eq. ([74|) by performing the 
change (a, b) (6, a) and noting that p{r) has even par- 
ity, i.e., p{— K,) = p{Hi), while p'{r) has odd parity, i.e., 
p\-k) = -p\i^). 
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IV. OTHER APPLICATIONS OF ELLIPTIC 
FUNCTIONS IN PHYSICS 

There are many more applications of elliptic func- 
tions in physics, including particle orbits and light 
paths in general relativity^i and solutions of cosmolog- 
ical models. These applications unfortunately involve 
advanced topics that fall well outside the purpose of the 
present work. 

One interesting application worthy of discussion, how- 
ever, deals with exact solutions of the Korteweg-de Vries 
(KdV) equationi^ 



du 
'dt 



du 
dx 







(75) 



which describes the nonlinear evolution of the field 
u(x^ t). This nonlinear equation appears in many areas of 
physic^ and is a member of an important class of non- 
linear partial differential equations that possesses soliton 
solutions.^^'^^'^^ 

A travelling-wave solution of the Korteweg-de Vries 
(KdV) equation ([75]) is a function of the form u{x^ t) = 
'u(^), where ^ = tz {x — ct) denotes the wave phase (with 
constants and c to be determined). Substituting this 
travelling-solution in the KdV equation, we obtain an 
ordinary differential equation for v{£^): 



{v - c) v' ^ v" 



0, 



which can be integrated with respect to ^ to yield 

1 o 



2 // I 

n V = a ^ cv 



(76) 



where a is a constant of integration. If we multiply 
Eq. ([76|) with v' and integrate again with respect to ^, 
we obtain 

^ {v'f = {av + 13) + lv' - \v\ (77) 

where /3 is a second constant of integration. It is now im- 
mediately clear that i;(^) = ^4 p(^) + 5 can be expressed 
in terms of elliptic functions (where A = —12 h?^ B = c) 
because the right side of Eq. ([77|) involves a cubic poly- 
nomial in V. The travelling- wave solution of the KdV 
equation ([75|) is 

u{x^t) = Ap[hz{x — ct) + 7] + 5, 

where the constant 7 is determined from the initial con- 
dition u{x^O) = uo{x). 

Using the relation between the Weierstrass and Jacobi 
elliptic functions (see Appendix [X]) , the travelling- wave 
solution to the KdV equation ([75|) may also be expressed 



as 
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u{x^t) = a cn^[/^ (x — ct) + 7 I m] + 6, 



where m = \/(r3 — r2 ) / {r^ — ^1 ) , a = ra — r2 , 6 = r2 , and 
Hi = \/ljs^^V^J/6; here, ri < r2 < rs are the roots of 



the cubic polynomial on the right side of Eq. ([77|) . This 
second representation is known as the periodic cnoidal- 
wave solution of the KdV equation ([75|) . 

Lastly, we note that for the special case a = = /3 in 
Eq. ([77|) , for which rs = 3 c and ri = = r2 , then we find 
m = l,a = 3c, 6 = 0,/^ = and the travelling- wave 

solution becomes 



u{x^t) = Scsech 



{x — ct) 



which describes the well-known localized soliton solution 
of the KdV equation ([75|) . 



V. SUMMARY 

We presented a brief introduction of the Jacobi and 
Weierstrass elliptic functions with applications in classi- 
cal mechanics. The problem of the planar pendulum was 
used to establish a connection between the Jacobi and 
Weierstrass elliptic functions. The double periodicity of 
the elliptic functions was easily observed in each of the 
problems of classical mechanics discussed in Sees. lU and 
Hm We also briefly discussed applications of the elliptic 
functions in other areas of physics (e.g., travelling- wave 
solutions of the Korteweg-de Vries equation). 

Because of their similarity with trigonometric func- 
tions sinz and cosz, physicists are most familiar with 
the Jacobi elliptic functions sn(2:|m) and cn(z|m). This 
familiarity is further increased by the fact that the prob- 
lems of the planar pendulum and the force- free asymmet- 
ric top are solved simply in terms of the Jacobi elliptic 
functions (as was demonstrated in Sec. [TT|). 

Physicists are usually less familiar with the Weierstrass 
elliptic function p{z; gs), introduced in Sec. IIIIl The 
complexity of the Weierstrass elliptic function is partly 
due to the fact that the form of its fundamental period 
parallelogram depends on the invariants ^2 and gs (and 
the modular discriminant A = ~ 27^1). This is in 
contrast to the fundamental period parallelogram of the 
Jacobi elliptic function which remains rectangular for all 
values of m 7^ 1. The greatest advantage of the Weier- 
strass elliptic function, however, is that the function p{z) 
itself and its derivative p^z) can be used to represent any 
doubly-periodic function f{z) = A{p) p' + B{p)^ where 
A and B are arbitrary functions. For example, using 
Eq. ([59]), we find sri^ {f^z\m) = (62 - 63)"^ [p{z;g2,g3) - 
63]. This simplicity was further demonstrated by the 
solution of the planar pendulum in terms of a single ex- 
pression z{r) = 1 — cos 0{r) = 2 p{r + 7) + /i for all values 
of energy. 

While the introduction of the Jacobi and Weierstrass 
elliptic functions in the undergraduate curriculum re- 
mains difficult, it is hoped that the present primer, com- 
bined with the easy access to mathematical software, will 
facilitate their dissemination. 



16 



APPENDIX A: RELATION BETWEEN 
WEIERSTRASS AND JACOBI ELLIPTIC 
FUNCTIONS 

In this Appendix, we explore the connection between 
the Jacobi and Weierstrass elhptic functions. For this 
purpose, we begin with the differential equation for the 
Weierstrass elliptic function 



dy . ^ . 

^ ) = 4y - 92y - g?>, 

and introduce the transformation 

y(^x) = a s^{z) + 



(Al) 



(A2) 



where p is an integer (7^ 0, 1), the function s{z) depends 
on the new coordinate z = nx^ and {a^jS^n) are con- 
stants to be determined. Under the transformation (|A2p . 
Eq. (|Aip becomes 



4a 



12/3 2 
s 



(12/32 



92 



+ (4/33 - ^2/3 - ^3) ^y^. (A3) 

The constants (a, /3, and the integer p ^ 0, 1 are chosen 
such that the right side of Eq. ()A3p has the Jacobi form 
(1 — s^) {1 — ms^). For p = 2, we recover the Jacobi form 
if 



(m + 1) 



(A4) 



is a root of the cubic polynomial 4/3^ — 92 P — gs (i.e., 

/3 = ei, 62, or 63) and 

2 12/32-ff2 
a = m K = . 

For p = —2, on the other hand, we recover the Jacobi 
form if /3, given by Eq. (|A4p . is again a root of the cubic 
polynomial 4/3^—^2/3 — ^3 (i-e., /3 = ei, 62, or 63) and 



a 



12/?"-ff2 
4 m /t:^ 



Hence, the Jacobi elliptic function s{z) is related to the 
Weierstrass elliptic function y{x) in the case of p = ± 2. 

An application of the first transformation {p = 2) 
shows that for m = (62 — 63) / (ei — 63) and = V^eT^^^, 
we find a = 62 — 63 and /3 = 63, and we obtain the relation 

p(x + cj2;^2,^3) = 63 + (62 - 63) cn^ (/^x I m) , (A5) 



or 



p(x + cc;3;^2,^3) = 63 + (62 - 63) sn^ (/^x I m) , (A6) 



which oscillate between 62 = p{oo2) = p{u;-\-u;s) and 63 = 
p{lj-\-u;2) = pi^s) (see Fig.[T6|). where a; = K{m)/i<i. Re- 
lation (|A6p plays a crucial role in expressing the Weier- 
trassian solution of the planar pendulum in terms of its 
Jacobian solution [see Eq. ([59])]. 

An application of the second transformation (jp = —2) 
shows that for the same m and as the first transforma- 
tion, we find Of = 61 — 63 and /3 = 63, and we obtain the 
relation 



p{x;g2,g3) 



63 



ei - 63 
sn^ {f<ix I m) ' 



(A7) 



which has singularities at x = and 2 K{m) / n = 2ix) and 
a minimum (61) at x = K{m)/n = co (see upper plot in 
Fig. [T3|). We note that the relation ()A7p is equivalent to 
the property (|48|) of the Weierstrass elliptic function: 



p{x) 



63 



(61 - 63) (62 - 63) 

p{x + iJs) - 63 



when substituting Eq. (|A6p into the denominator of 
Eq. (|A7l) . 



APPENDIX B: MATHEMATICAL DETAILS 

In this Appendix, we present mathematical details that 
relate the half-periods (|46|) of the Weierstrass elliptic 
function with the quarter periods K and K' of the Jacobi 
elliptic function. We assume throughout this Appendix 
that ^3 > and A > (i.e., < e < 1). 

First, we start with cc;i(^2,^3), where 5 > 61 > 62 > 
63 in the integral (|46|) . By introducing the change of 
variable^^ 5 = 63 + (61 — 63) csc^V^, we readily obtain 



^1(^2,^3 



n7r/2 



(61 - 63) - (62 - 63) sin^ i/j 



K{m) 

V^l - 63 



^(^2,^3), 



(Bl) 



where u is defined in Eq. (|42|). the modulus m{(j)) of the 
Jacobian quarter period is defined by the relation 



m(0) 



62 - 63 _ sin[(7r - 0)/3] 
61 - 63 sin[(7r + 0)/3] ' 



using the definitions dH]), and < ^ < tt along the seg- 
ment (b) in Fig. [T2l 

Next, we look at 6^3(^2, ^'3), which can also be expressed 

as 



^3(^2,^3) 



ds 



- f 

2 i-c 



i/4s3 - g^s- g3 
ds 



(B2) 



A/(ei - s)(e2 - s)(e3 - s) 



17 



By introducing the change of variable s = ei — {ei 
63) csc^?/^, we readily obtain 



^3(^2,^3) = i 



7t/2 



^J (ei - 63) - (ei - 62) sin^ i/j 
iK{m') 



^'(^2,^3), 



V^l - 63 

where cj' is defined in Eq. (|43|) and 

m'(0) = 1 — m(^) = ^"^ 



(B3) 



sin(0/3) 



ei - 63 sin[(7r + (/))/3] * 
Lastly, we look at ^^2(^2,^3), which can be written as 



002 = ^1 



^ 62 



(is 



^"(61 - s){s - e2){s - 63) 



(B4) 



We now introduce the change of variable^ 



^1 _ _ (ei - 62) (62 - 63) 



(s - 62) 



and readily find that the second integral in Eq. (|B4j) is 
simply uj 3 (92,93) so that 



^2(^2,^3) = ^1(^2,^3) + ^3(^2,^3) 



a; + a; 



Vei - 63 * 



(B5) 
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